% Hill_elastance.m
% function calculating time-varying elastance curve
% approximation as described by Stergiopulos et al., Am J Physiol 1996
% double Hill function

function [fHill] = elastance(tN,amp_Hill,t01_Hill,t02_Hill,n_Hill,m_Hill)

% default parameter values

fHill=amp_Hill*(((tN/t01_Hill).^n_Hill./(1+(tN/t01_Hill).^n_Hill)).*(1./(1+(tN/t02_Hill).^m_Hill)));